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Abstract. Let A* > denote the largest possible value of A such that 
A^u = Ae" in B 

du 

u = — = on oB 

dn 

has a solution, where B is the unit ball in and n is the exterior unit normal 
vector. We show that for A = A* this problem possesses a unique weak solution 
u* . We prove that u* is smooth if Af < 12 and singular when N > 13, in which 
case u*{r) = -41ogr + log(8(Af - 2)(Af - 4)/A*) + o(l) as r 0. We also 
consider the problem with general constant Dirichlet boundary conditions. 



1. Introduction 
We study the fourth order problem 



(1) 



A^u = Ae" in B 
u = a on dB 

^— ~h on dB 



dn 

where a, 6 G R, B is the unit baU in K", iV > 1, n is the exterior unit normal 
vector and A > is a parameter. 

Recently higher order equations have attracted the interest of many researchers. 
In particular fourth order equations with an exponential non-linearity have been 
studied in 4 dimensions, in a setting analogous to Liouville's equation, in [3, 12, 24] 
and in higher dimensions by [1, 2, 4, 5, 13]. 

We shall pay special attention to (1) in the case a = 6 = 0, as it is the natural 
fourth order analogue of the classical Gelfand problem 

Aw = Ae" in O 
u = on dVl 



(2) 



(n is a smooth bounded domain in R^) for which a vast literature exists [7, 8, 9, 
10, 18, 19, 20, 21]. 

From the technical point of view, one of the basic tools in the analysis of (2) is 
the maximum principle. As pointed out in [2], in general domains the maximum 
principle for A^ with Dirichlet boundary condition is not valid anymore. One of the 
reasons to study (1) in a ball is that a maximum principle holds in this situation, 
see [6] . In this simpler setting, though there are some similarities between the two 
problems, several tools that are well suited for (2) no longer seem to work for (1). 
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As a start, let us introduce the class of weak solutions we shall be working with : 
we say that u G H^{B) is a weak solution to (1) if e" G L^{B), u = a on dB, 



§^ = hondB and 

on 

[ AuAip = x[ e"ip, for aU </3 G C^(B). 
Jb Jb 

The following basic result is a straightforward adaptation of Theorem 3 in [2] . 

Theorem 1.1. ([2]) There exists X* such that i/O < A < A* then (1) has a minimal 
smooth solution ux and if X > X* then (1) has no weak solution. 

The limit u* = lim;^^;^» ux exists pointwise, belongs to H^{B) and is a weak 
solution to (1). It is called the extremal solution. 

The functions ux, < X < X* and u* are radially symmetric and radially de- 
creasing. 

The branch of minimal solutions of (1) has an important property, namely ux is 
stable in the sense that 



(3) / (A^)^> A / e"V^ V^GCo°°(i?), 

J B J B 

see [2, Proposition 37]. 

The authors in [2] pose several questions, some of which we address in this work. 
First we show that the extremal solution u* is the unique solution to (1) in the class 
of weak solutions. Actually the statement is stronger, asserting that for A = A* 
there are no strict super-solutions. 

Theorem 1.2. // 

(4) V G H\B), G L\B), v\ob = a, f^las < b 
and 

(5) / AvAip>X* [ e> V</7GCo°°(B), (/?>0, 
Jb Jb 

then V ~ u* . In particular for A = A* problem (1) has a unique weak solution. 

This result is analogous to work of Martel [19] for more general versions of (2) 
where the exponential function is replaced by a positive, increasing, convex and 
superlinear function. 

Next, we discuss the regularity of the extremal solution u* . In dimensions = 
5, . . . , 16 the authors of [2] find, with a computer assisted proof, a radial singular 
solution Ua to (1) with a = 6 = associated to a parameter Xa > 8(A — 2)(A^ — 4). 
They show that Ao- < A* if A^ < 10 and claim to have numerical evidence that 
this holds for A < 12. They leave open the question of whether u* is singular in 
dimension A < 12. Wc prove 

Theorem 1.3. If N < 12 then the extremal solution u* of (1) is smooth. 

The method introduced in [10, 20] to prove the boundedness of u* in low di- 
mensions for (2) seems not useful for (1), thus requiring a new strategy. A first 
indication that the borderline dimension for the boundedness of u* is 12 is Rellicli's 
inequality [23] , which states that if A^ > 5 then 



16 yR^R 



(6) (A^r> \' V^gCo-(R^), 



STABLE SOLUTIONS FOR A^ti = Ae 



3 



where the constant N'^{N — 4)^/16 is known to be optimal. The proof of Theo- 
rem 1.3 is based on the observation that if u* is singular then A*e" ~ 8{N —2){N — 
4)|a;|-4 near the origin. But 8(A^ - 2){N - 4) > N^{N - 4)Vl6 if iV < 12 which 
would contradict the stability condition (3). 

In view of Theorem 1.3, it is natural to ask whether u* is singular in dimension 

> 13. If a = 5 = 0, we prove 

Theorem 1.4. Let N > 13 and a = b = 0. Then the extremal solution u* to (1) 
is unbounded. 

For general boundary values, it seems more difficult to determine the dimensions 
for which the extremal solution is singular. We observe first that given any a, 6 S R, 
u* is the extremal solution of (1) if and only if w* — a is the extremal solution of 
the same equation with boundary condition m = on dB. In particular, if A* (a, b) 
denotes the extremal parameter for problem (1), one has that X*{a,b) — e~'^A*(0,6). 
So the value of a is irrelevant. But one may ask if Theorem 1.4 still holds for any 
TV > 13 and any 6 S R. The situation turns out to be somewhat more complicated : 

Proposition 1.5. 

a) Fix iV > 13 and take any a e M. Assume b > —4. There exists a critical 
parameter 5"""^ > 0, depending only on N , such that the extremal solution 
u* is singular if and only if b < 5"'°^. 

b) Fix b > —4 and take any a £ R. There exists a critical dimension iV™*" > 
13, depending only on b, such that the extremal solution u* to (1) is singular 
if N> iV™". 

Remark 1.6. 

• We have not investigated the case b < —4. 

• If follows from item a) that for b £ [—4, 0], the extremal solution is singular 
if and only if N > 13. 

• It also follows from item a) that there exist values ofb for which N™™ > 13. 
We do not know whether u* remains bounded for li < N < iV™™. 

Our proof of Theorem 1.4 is related to an idea that Brezis and Vazquez ap- 
plied for the Gelfand problem and is based on a characterization of singular energy 
solutions through linearized stability (see Theorem 3.1 in [8]). In our context we 
show 

Proposition 1.7. Assume that u e H^{B) is an unbounded weak solution of (1) 
satisfying the stability condition 

(7) A / < f (^^)2^ ^ ^oo^j^y 

J B J B 

Then A = A* and u = u* . 

Wc do not use Proposition 1.7 directly but some variants of it - see Lemma 2.6 
and Remark 2.7 in Section 2 - because we do not have at our disposal an explicit 
solution to the equation (1). Instead, we show that it is enough to find a sufficiently 
good approximation to u* . When TV > 32 we are able to construct such an ap- 
proximation by hand. However, for 13 < TV < 31 we resort to a computer assisted 
generation and verification. 
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Only in very few situations one may take advantage of Proposition 1.7 directly. 
For instance for problem (1) with a = and b = —4 we have an explicit solution 



associated to A = 8{N — 2){N — 4). Thanks to Rellich's inequality (6) the solution 
u satisfies condition (7) when N > 13. Therefore, by Theorem 1.3 and a direct 
application of Proposition 1.7 we obtain Theorem 1.4 in the case b = —4. 

In [2] the authors say that a radial weak solution u to (1) is weakly singular if 



For example, the singular solutions of [2] verify this condition. 
As a corollary of Theorem 1.2 we show 

Proposition 1.8. The extremal solution u* to (1) with b > —4 is always weakly 
singular. 

A weakly singular solution cither is smooth or exhibits a log-type singularity at 
the origin. More precisely, if it is a non-smooth weakly singular solution of (1) with 
parameter A then (see [2]) 



In Section 2 we describe the comparison principles we use later on. Section 3 is 
devoted to the proof of the uniqueness of u* and Propositions 1.7 and 1.8. We prove 
Theorem 1.3, the boundcdness of u* in low dimensions, in Section 4. The argument 
for Theorem 1.4 is contained in Section 5 for the case > 32 and Section 6 for 
13 < A^ < 31. In Section 7 we give the proof of Proposition 1.5. 

Notation. 

• Br: ball of radius R in centered at the origin. B = Bi. 

• n: exterior unit normal vector to Bfi 

• All inequalities or equalities for functions in spaces are understood to 



u{x) 



41og|x| 



Vim. ru'(r) exists. 



lim u{r) -t- 4 log r = log 



%{N -2){N -4) 
A 



lim ru'{r) = —4. 



be a.e. 



2. Comparison principles 



Lemma 2.1. (Boggio's principle, [6]j If u G C'^{Bji) satisfies 

( A^u > in Br 




u = — =0 on dBn 



then u > in Br . 
Lemma 2.2. Let u E 



{Bji) and suppose that 




for all ip £ C^{Bii) such that ip > in Br, (p\dB 
Br. Moreover u = or u > a.e. in Br. 



R 



= ^\dBR- Then u > in 



For a proof see Lemma 17 in [2]. 
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Lemma 2.3. If u ^ H^[Br) is radial, A^u > m Br in the weak sense, that is 

I AuA<^>0 G C^(Sfl), > 
Jbr 

and ulasH > 0, f^|9B„ < then u > in Br. 

Proof. Wc only deal with the case i? = 1 for simplicity. Solve 

A^wi = A^M in Bi 

-Ul = — — = on oBi 
on 

in the sense ui G H^{Bi) and AmiA(p = AuA(/9 for all G C^{Bi). Then 
Ml > in i?i by Lemma 2.2. 

Let U2 = u — Ul so that A^U2 = in Bi. Define / = Au2. Then A/ = in Bi 
and since / is radial we find that / is constant. It follows that U2 — ar^ + b. Using 
the boundary conditions we deduce a + b > and a < 0, which imply U2 > 0. □ 

Similarly we have 

Lemma 2.4. If u G H^{Br) and A^u > in Bp in the weak sense, that is 
[ AuAif > G C^{Br), ip>Q 

J Br 

and u\oBa = 0, ^{oBr < then u > in Br. 

The next lemma is a consequence of a decomposition lemma of Moreau [22] . For 
a proof see [14, 15]. 

Lemma 2.5. Let u G Hq{Br). Then there exist unique w,v G H^{Br) such that 
u ~ w + V, w > 0, A^u < in Br and J^^ AwAv = 0. 

We need the following comparison principle. 

Lemma 2.6. Let ui, U2 G H^{Br) with e"\ e"^ G L^{Br). Assume that 

A^ui < Ae"i m Br 

in the sense 

(8) f AuiAip<xf e"V ytpeC^\BR),ip>0, 

J Br J Br 

and A^U2 > Ae"^ in Br in the similar weak sense. Suppose also 

, dui du2 1 

ui\dBR = U2\dBR and ^— as^ = ^r— asH- 

On an 

Assume furthermore that ui is stable in the sense that 

(9) a / e"V'</ (A^)^ ^^eC^{BR). 



Then 

Ul < U2 in Br. 
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Proof.Let u = ui — U2- By Lemma 2.5 there exist w,v E Hq{Bii) such that 
u = w + v,w>0 and A^w < 0. Observe that v < so w > ui — U2- 
By hypothesis wc have for all ip £ C^{Bji), ip> 0, 



Br J Br, J BRn[ui>ui] 

and by density this holds also for w. 
(10) / [Awf ^ I A{ui~U2)Aw 

Jbr JBr 



J BRn[ui>U2] Jbr 

where the first equality holds because AwAv = 0. By density we deduce from 
(9): 

(11) A j e"iw2 < [ {Awf. 

JBr J Br 

Combining (10) and (11) we obtain 

Br J Br 

Since ui — U2 < w the previous inequality implies 

(12) < / (e"i - - e"i {ui - U2))w. 

JBr 

But by convexity of the exponential function e"^ — e""^ — e"i(ui — U2) < and we 
deduce from (12) that (e"^ — e"^ — e"i {ui — U2))w ~ 0. Recalling that ui — U2 < w 
we deduce that ui < U2. □ 

Remark 2.7. The following variant of Lemma 2.6 also holds: 

Let ui, U2 e H^{Bii) be radial with e"i, e^^ G L'^{Bb)- Assume A'^ui < Ae^i 
in Bpt in the sense of (8) and A'^U2 > Ae"^ m B/j. Suppose ui\qbr < U2|a_Bp anc? 
^^laSn — ^^I^Sh aJ^rf ^/ifl^ the stability condition (9) holds. Then ui < 112 in Br. 

Proof. We solve for u e H'^{Bji) such that 

/ AuAip = / A(ui - u2)Aip \fip e C^{Br). 

JBr J Br 

By Lemma 2.3 it follows that u > ui — U2. Next we apply the decomposition of 
Lemma 2.5 to u, that isu^w + v with w, tj e H^{Bb). ui > 0, A^v < in Bji and 
AtiiAu = 0. Then the argument follows that of Lemma 2.6. □ 

Finally, in several places we will need the method of sub and supersolutions in 
the context of weak solutions. 

Lemma 2.8. Let A > and assume that there exists u G H^{Bfj) such that 
e° e L\Br), 

I AuAip >\ I e"<p for all ip e C^{BR),ip > 
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du 

u = a, — < on oBi. 
an 

Then there exists a weak solution to (1) such that u < u. 
The proof is similar to that of Lemma 19 in [2]. 

3. Uniqueness of the extremal solution: proof of Theorem 1.2 

Proof of Theorem 1.2. Suppose that v e H'^{B) satisfies (4), (5) and v ^ u* 
Notice that we do not need v to be radial. 
The idea of the proof is as follows : 



Step 1. The function 



"0 = + v) 



2 

is a super-solution to the following problem 

A^u = A*e" + ^ir|e^ in B 



(13) 



u = a 



on dB 



^— = b on dB 

dn 



for some /i = /iq > 0, where rj E C^{B), < 77 < 1 is a fixed radial cut-off function 
such that 

r]{x) = 1 for \x\ < r]{x) = for |a;| > |. 

Step 2. Using a solution to (13) we construct, for some A > A*, a super-solution 
to (1). This provides a solution u\ for some A > A*, which is a contradiction. 

Proof of Step 1. Observe that given < i? < 1 we must have for some cq = 
co(i?) > 

(14) v{x) >u*{x) +ca \x\<R. 

To prove this we recall the Green's function for with Dirichlet boundary condi- 
tions 

AlG{x, y) = Sy xeB 
G{x, y)=0 xedB 
dC 

— {x,y)=0 xedB, 

where Sy is the Dirac mass at y e _B. Boggio gave an explicit formula for G(x, y) 
which was used in [16] to prove that in dimension N > 5 (the case 1 < < 4 can 
be treated similarly) 

d(x)2d(y)2' 



(15) Gix,y)^\x-y\^-^min[l 



\x-y\ 

where 

d{x) = dist{x, OB) ^l~\x\ 



4 
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and a ^ b means that for some constant C > we have C < b < Ca (uniformly 
for x,y E B). Formula (15) yields 

(16) G{x,y)>cd{x)^d{yf 

for some c > and this in turn implies that for smooth functions v and u such that 
v-ne H^iB) and A'^{i - u) > 0, 

-/ \ \ f fdA^G, . ^, .d{v — u)\ , 

v{y) - u{y) = J \^-^—{x,y)[v - u) - A^G(x,y) — — — j dx 

+ / G{x,y)A'^{v -u)dx 
Jb 

> cd{yf I [A^v - A^u)d[xf dx. 
Jb 

Using a standard approximation procedure, we conclude that 

v{y)-u*{y)>cd{yfX* [ {e" - e""' )d{xf dx. 



Since v > u* , v ^ u* we deduce (14). 

Let uq = {u* + v)/2. Then by Taylor's theorem 

(17) e" = e"" + {v- uo)e"» + - uo)'e"« + - uo fe^" + ^iv- «o)^e«^ 

2 O z4 

for some Uq < S,2 l£ v and 
(18) 

e"* = e"" + {u* - uo)e"" + l:{u* - uo)'e"" + W ~ "o)'e"° + ^{u* - uo fe^' 

2 z4 

for some m* < < uq. Adding (17) and (18) yields 

(19) i(e^+e"')>e"° + -(w-u*)^e'"'. 

2 8 

From (14) with R = 3/4 and (19) we see that uq = {u* + v)/2 is a super-solution 

of (13) with /iQ := co/8. 

Proof of Step 2. Let us show now how to obtain a weak super-solution of (1) for 
some A > A*. Given fi > 0, let u denote the minimal solution to (13). Define (pi as 
the solution to 

A'^ipi = fj.rje'^ in B 
(fi = on dB 

dip 



and ip2 be the solution of 



^ i = on dB, 
on 



AV2 =0 in B 

Lp2 = a on dB 

-- — = on dB. 

dn 



If > 5 (the case 1 < iV < 4 can be treated similarly), relation (16) yields 
(20) ipi{x) > cid(x)^ for aU x G B, 
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for some ci > 0. But it is a radial solution of (13) and therefore it is smooth in 
B\Bi/i. Thus 

(21) u{x) < Mipi + ip2 for all x € B1/2, 

for some AI > 0. Therefore, from (20) and (21), for A > A* with A - A* sufficiently 
small we have 

(^ - l)u < + (^ - 1)(^2 inS. 

Let w = -^u — ipi — (p- — l)(p2- The inequality just stated guarantees that w < u. 
Moreover 

.2„.. , V, 



and 



A^w = Ae" + — r;e" - firje" > Ae" > Ae"' in B 
A* 



w = a ~ b on dB. 

on 



Therefore w is & super-solution to (1) for A. By the method of sub and super- 
solutions a solution to (1) exists for some A > A*, which is a contradiction. □ 

Proof of Proposition 1.7. Let u £ H^{B), A > be a weak unbounded solution 
of (1). If A < A* from Lemma 2.6 we find that u < u\ where u\ is the minimal 
solution. This is impossible because u\ is smooth and u unbounded. If A = A* then 
necessarily u = u* by Theorem 1.2. □ 

Proof of Proposition 1.8. Let u denote the extremal solution of (1) with h > —4. 
If u is smooth, then the result is trivial. So we restrict to the case where u is 
singular. By Theorem 1.3 we have in particular that N > 13. We may also assume 
that a = 0. If 6 = -4 by Theorem 1.2 we know that if TV > 13 then u = -4 log \x\ 
so that the desired conclusion holds. Henceforth we assume b > —4 in this section. 
For p > define 



so that 



Then 



Up{r) = u{pr) + 41og/9, 
A^Mp = A*e"'- in B^/p. 



dp 

Hence, there is J > such that 



= u'(l) + 4>0. 

p— l,r— 1 



Up{r) < u{r) for all 1 - (5 < r < 1, 1 - (5 < p < 1. 

This implies 

(22) Up{r)<u{r) for all < r < 1, 1 - (5 < p < 1. 
Otherwise set 

7'o = sup { < 7' < 1 1 Up{r) > u(t) }. 

This definition yields 

(23) Mp(ro)=u(ro) and Up(?'o) < u'(r'o). 
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Write a = u(ro), (3 = u'{ro). Then u satisfies 







= Ae 


(24) 


u{ro) - 


= a 




. "'(''o) = 


= 13. 



(r) = u'{r)r + 4 for all < r < 1 

p=i 



on Brg 



Observe that u is an unbounded H^{Brg) solution to (24), which is also stable. 
Thus Proposition 1.7 shows that u is the extremal solution to this problem. On the 
other hand Up is a supersolution to (24), since Up(?'o) £ /9 by (23). We may now 
use Theorem 1.2 and we deduce that 

u{r) = Up{r) for all < r < Tq, 

which in turn implies by standard ODE theory that 

u{r) ~ Up{r) for all < r < 1, 

a contradiction with (22). This proves estimate (22). 
From (22) we see that 

nil 

(25) — ^ (r) > for all < r < 1. 

dp p=i 

But 

dup 
dp 

and this together with (25) implies 

(26) ^(r) = -{u'{pr)pr + 4) > for all < r < -, < p < 1. 

dp p p 

which means that itp(r) is non-decreasing in p. We wish to show that limp^o Up{r) 
exists for all < r < 1. For this we shall show 

(27) Up{r) > -41og(r) + log ^ ^^^ ~ ^^J^ ~ for all < r < i, < p < 1. 
Set 

' 8{N - 2){N - iy 
A^^ 

and suppose that (27) is not true for some < p < 1. Let 

ri = sup { < r < 1/p I Up{r) < uo{r) }. 

Observe that 

(28) A* > 8(7V-2)(iV-4). 

Otherwise w = — 41nr would be a strict supersolution of the equation satisfied by 
u, which is not possible by Theorem 1.2. In particular, ri < 1/p and 

Up{ri) = ito(r-i) and u'p{ri) > -Uo(ri). 

It follows that Mo is a supersolution of 



■"o('') = -41og(r) + log 



(29) 





= A*e 


< ^ 




du 






= B 


^ dn 



in 

on dBr^ 



on dBr^ , 
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with A ~ Up{ri) and B = u'p{ri). Since Up is a singular stable solution of (29), it 
is the extremal solution of the problem by Proposition 1.7. By Theorem 1.2, there 
is no strict supersolution of (29) and we conclude that Up = uq first for < r < ri 
and then for < r < 1/p. This is impossible for p > because Up{l/p) ~ 41ogp 
and uo(l/p) < 41ogp + log( ^^^~^^l^^'^^ ) < Up(l/p) by (28). This proves (27). 
By (26) and (27) we see that 

v{r) = lim Up{r) exists for all < r < +cx3, 

where the convergence is uniform (even in for any fc) on compact sets of M^\{0}. 
Moreover v satisfies 

(30) A^z; = X*e" in \ {0}. 
Then for any r > 

v{r) = lim Up{r) = lim u{pr) + 41og(pr) — 41og(r) = v{\) — 41og(r). 

Hence, using equation (30) we obtain 

, /8(iV-2)(7V-4)\ 
v{r) = -41ogr + log ( ^ \ = uo(0- 

But then 

u'p(r) — u'{pr)p —4, as p 0. 
and therefore, with r = 1 

(31) pu'{p)^-A asp^O. 

□ 



4. Proof of Theorem 1.3 

We will show first 

Lemma 4.1. Suppose that the extremal solution u* to (1) is singular. Then for 
any tr > there exists < i? < 1 such that 

(32) u*{x)>{l-a)\og(^-^y \f\x\<R. 

Proof. Assume by contradiction that (32) is false. Then there exists cr > and a 
sequence Xk G B with — !■ such that 

(33) u*{xk)<{l-cT)\og(^ ^ 
Let Sfc = \xk\ and choose < Afc < A* such that 



(34) max = (0) = log ( ^ 

B \4. 

Note that — > A*, otherwise would remain bounded. Let 

= - — -y— X e Bk ^ —B. 
log -V) Sk 
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Then 0<vk<l, Vk{0) = 1, 

log pr 

by (34). By elliptic regularity Vk v uniformly on compact sets of to a 
function v satisfying < u < 1, w(0) = 1, A^w = in R^.By Liouville's theorem 
for biharmonic functions [17] we conclude that v is constant and therefore v = 1. 
Since l^fej = Sfe we deduce that 

log(i) ' 

k 

which contradicts (33). □ 

Proof of Theorem 1.3. We write for simplicity u = u* , A = A*. Assume by 
contradiction that u* is unbounded and 5 < < 12. If < 4 the problem is 
subcritical, and the boundedness of u* can be proved by other means : no singular 
solutions exist for positive A (see [2])-though in dimension A^ = 4 they can blow up 
as A — > 0, see [24]. 

For e > let ip = and let t] e C^(R^) with = 1 in and 

supp{rj) C B. Observe that 

{Ai'f = {Hn + 0(e))|x|-^+2e^ ^^^^^ ^ ^ ) 

Using a standard approximation argument as in the proof of Lemma 2.6, we can 
use ipr] as a test function in (9) and we obtain 



(A^/.)' + 0(1) > A / e"V2^ 



B 

since the contribution of the integrals outside a fixed ball around the origin remains 
bounded as e — ^ (here 0(1) denotes a bounded function as e ^ 0). 
This implies 

(35) A / e"|x|4-™ < {Hn + 0(e)) / = uj^^ + 0(1) 

Jb Jb ^£ 

where ujn is the surface area of the unit A — 1 dimensional sphere S*^^^. In 
particular Jq e"\x\^~^^'^^ < +00. 

For e > we define (p — Note that away from the origin 

(36) AV = efc^r|x|~^+2^ where fc^v = 4(A - 2)(Af - 4) + 0(e). 
Let ipj solve 

r A^ip^ = £fcJvmin(|x|-^+2^ j) in B 

(37) <^ d^, 
\ ^j = ^^0 on dB. 



d 



n 
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Then ipj '\ Lp as j ^ +00. Using (35) and (37) 

ekN I u min(|x|-^+2^j) = ( uAVj = A / e^ipj 

JB JB JB 



< A / e"^ 



B 



2e 

where 0(1) is bounded as e independently of j. Letting j — > +00 yields 

(38) ekN I u < iON^ + 0(1), 

Jb 2e 

showing that the integral on the left hand side is finite. On the other hand, by (32) 

(39) ekN I u |a;r^+2e > skNUJN{l - cr) f \og{^)r-^+^' dr = fcjvwjv(l - a)-. 

JB Jo T ^ 

Combining (38) and (39) we obtain 

{l~a)kN<^ + 0{e). 



Letting e ^ and then cr ^ we have 

8{N^2){N~A)<Hn 

16 

This is valid only if > 13, a contradiction. □ 

Remark 4.2. The conclusion of Theorem 1.3 can he obtained also from Proposi- 
tion 1.8. However that proposition depends crucially on the radial symmetry of the 
solutions, while the argument in this section can be generalized to other domains. 

5. The extremal solution is singular in large dimensions 

In this section we take a = 6 = and prove Theorem 1.4 for N > 32. 

The idea for the proof of Theorem 1.4 is to to estimate accurately from above the 
function A*e" , and to deduce that the operator — A*e'" has a strictly positive 
first eigenvalue (in the Hq{B) sense). Then, necessarily, u* is singular. 

Upper bounds for both A* and u* are obtained by finding suitable sub and 
supersolutions. For example, if for some Ai there exists a supersolution then A* > 
Ai. If for some A2 one can exhibit a stable singular subsolution u, then A* < A2. 
Otherwise A2 < A* and one can then prove that the minimal solution ux^ is above 
u, which is impossible. The bound for u* also requires a stable singular subsolution. 

It turns out that in dimension iV > 32 we can construct the necessary subso- 
lutions and verify their stability by hand. For dimensions 13 < iV < 31 it seems 
difficult to find these subsolutions explicitly. We adopt then an approach that 
involves a computer assisted construction of subsolutions and verification of the 
desired inequalities. We present this part in the next section. 

Lemma 5.1. Assume N > 13. Then u* < u ~ — 41og|x| in Bi. 
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(40) 



Proof. Define u{x) = — 41og|2:|. Then u satisfies 

A2?I = 8(iV-2)(iV-4)e" in M'^ 
u = on dBi 

^ = -4 on dBi 

on 

Observe tliat since w is a supersolution to (1) with a = 6 = we deduce imme- 
diately that A* > 8(A^ - 2)(iV - 4). 

In the case A* = 8{N — 2){N — 4) we have u\ < u for all < A < A* because u is 
a supersolution, and therefore u* < u holds. Alternatively, one can invoke Theorem 
3 in [2] to conclude that we always have A* > 8(A^ - 2)(A^ - 4). 

Suppose now that A* > 8(7V - 2)(A^ - 4). We prove that ux < u for all 8{N - 
2){N — 4) < A < A*. Fix such A and assume by contradiction that ux < u is not 
true. Note that for r < 1 and sufficiently close to 1 we have ux{r) < u{r) because 
u'y^{l) = while ft'(l) = -4. Let 

Ri = inf{ < i? < 1 I itA < M in (i?, 1) }. 

Then < i?i < 1, ux{Ri) = u{Ri) and u'-)^{Ri) < u'{Ri). So ux is a super-solution 
to the problem 

A^u = 8{N -2){N - 4)e" in Br, 
u^ux{Ri) ondBu, 
du 

— (i?i) ondBR, 

while u is a subsolution to (40). Moreover it is stable for this problem, since from 
Rellich's inequality (6) and 8(iV - 2)(A^ - 4) < N^{N -4)^/16 for N > 13, we have 

8(iV-2)(7V-4) / eV< / FR^/ (^^)' y^eC^iBR,). 

By Remark 2.7 we deduce u < in Br, which is impossible. □ 

An upper bound for A* is obtained by considering again a stable, singular sub- 
solution to the problem (with another parameter, though): 

Lemma 5.2. For N > 32 we have 

(41) A* < 8(A^-2)(7V-4)e^ 

Proof. Consider w = 2(1 — r^) and define 



where u{x) = — 41og|x|. Then 

A^u = 8(N - 2)(N - 4)4 = 8(N - 2)(N - 4)e" = 8(iV - 2)(N - 4)e"+"' 

<8{N ~2){N ~ 4:)e^e". 

Also u(l) = u'{l) = 0, so u is a subsolution to (1) with parameter Aq = 8{N — 
2)(Ar-4)e2. 

For > 32 we have Ao < N'^{N- 4)^/16. Then by (6) u is a stable subsolution 
of (1) with A = Aq. If a* > Ao = 8(7V-2)(A-4)e2 the minimal solution ux^ to (1) 
with parameter Aq exists and is smooth. From Lemma 2.6 we find u < uxg which 
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is impossible because u is singular and is bounded. Thus we have proved (41) 
for > 32. □ 

Proof of Theorem 1.4 in the case N > 32. 

Combining Lemma 5.1 and 5.2 we have that if TV > 32 then A*e"'" < r-'^ S{N - 
2) (TV - 4)e2 < r-^N^{N - 4)^/16. This and (6) show that 

which is not possible if u* is bounded. □ 



6. A COMPUTER ASSISTED PROOF FOR DIMENSIONS 13 < < 31 



Throughout this section we assume a = b = 0. As was mentioned in the previous 
section, the proof of Theorem 1.4 relies on precise estimates for u* and A*. We 
present first some conditions under which it is possible to find these estimates. 
Later we show how to meet such conditions with a computer assisted verification. 

The first lemma is analogous to Lemma 5.2. 

Lemma 6.1. Suppose there exist e > 0, A > and a radial function u G H^{B) fl 



loc \ {0}) such that 

A^t < Ae" 



Hl)\<e, 
u(^L^{B) 



for alio < r <1 
du , 



dn 



(1) 



< £ 



(42) 
Then 

Proof. Let 

(43) 
so that 

and 

It follows that 



Ae 



e Lp 



< 



{ALpf for all if e C^{B). 



A* < Ae 



2e 



ip{r) 



er 



2e 



A^V' = 0, V'(l) = ~£, V''(l) = 2e 
-2e < ip{r) < -e for aU < r < 1. 
A^{u + V') < Ae" = Ae-'^e"+'^ < Ae2^e"+'''. 



On the boundary we have u{l) + ipil) < 0, w'(l) + '0'(1) ^ 0- Thus u + ip is a 
singular subsolution to the equation with parameter Ae^'^. Moreover, since ip ^ ~£ 
we have Ae^^e"^''' < Ae^e" and hence, from (42) we see that u + ip is stable for the 
problem with parameter Ae^^. If Ae^^ < A* then the minimal solution associated to 
the parameter Ae^"^ would be above u + 'ip, which is impossible because u is singular. 
□ 
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Lemma 6.2. Suppose we can find e > 0, X>0 and ue H^iB) n Wi^^{B \ {0}) 
such that 

A^u > Ae" for all < r < 1 



1^(1) 
on 



< £. 



Then 



Xe 



-2e 



< X* 



Proof. Let ■0 be given by (43). Then u — -i/' is a supersolution to the problem with 
parameter Xe~^^. □ 

The next result is the main tool to guarantee that u* is singular. The proof, as 
in Lemma 5.1, is based on an upper estimate of u* by a stable singular subsolution. 

Lemma 6.3. Suppose there exist so,s > 0, Aa > and a radial function u G 
H^{B) n Wl^^^{B \ {0}) such that 



(44) 


A^U < {Xa 


+ 


£o)e" 


for all < 


r < 1 


(45) 


A^U > {Xa 




£o)e" 


for all < 


r < 1 


(46) 








< e 




(47) 


u^L°°{B) 










(48) 


Po f e V^ 
Jb 


< 


[ {Aipf for all 
Jb 


^ e Co"(i?), 


where 












(49) 




{Xa+eo) 

{Xa - eoY 




Then u* 


is singular and 










(50) 


{K- 


£0 


)e-^' <X* < {Xa + 





Proof. By Lemmas 6.1 and 6.2 we have (50). Let 

Aa + £0 



6 = log 



Aa — £0 



3e. 



and define 

ip{r) = -Ir" + 2(5. 
We claim that 

(51) u* <u + ip in Si. 
To prove this, we shall show that for A < A* 

(52) u\<u + (p in Bi. 
Indeed, we have 

AV -62N{N + 2) 

ip{r) > S for all < r < 1 

'p{l)>S>e, ip'{l)^-6< 
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and therefore 

A^u + ^)< {Xa + eo)e" + < (Aa + eo)e" = (A^ + eo)e-'^e"+'^ 

(53) < (Aa+£o)e-V+'^. 
By (50) and the choice of 5 

(54) (A,+£o)e-' = (Aa-£o)e-'^<A*. 

To prove (52) it suffices to consider A in the interval (Aa — £0)6""^^ < A < A*. Fix 
such A and assume that (52) is not true. Write 

u = u + if 

and let 

Ri = sup{ < i? < f I ux{R) = u{R) }. 

Then < i?i < f and ux{Ri) = u{Ri). Since u'^{l) = and u'{l) < we must 
have u'^{Ri) < u'{Ri). Then u\ is a solution to the problem 

A^u = Ae" in Br, 

u = ux{Ri) on dBji^ 

3u 

— = u'y^iRi) on SBr^ 
on 

while, thanks to (53) and (54), m is a subsolution to the same problem. Moreover 
u is stable thanks to (48) since, by Lemma 6.f , 

(55) A< A* < (Aa+£o)e'" 
and hence 

Ae" < (Aa+£o)e'^'e2*e" </3oe". 

We deduce m < ua in Bn^ which is impossible, since u is singular while ux is 
smooth. This establishes (51). 
From (51) and (55) we have 

A*e"* < /3oe-'^e" 

and therefore 

f„(A</j)2-A*e"V2 
inf , — > 0. 

This is not possible if u* is a smooth solution. □ 

For each dimension 13 < < 31 we construct u satisfying (44) to (48) of the 
form 



(56) u{r) = 



-41ogr + log r'^~Y'^-"M forO<r<ro 



A 

u(r) for ^0 < ?' < 1 

where ii is explicitly given. Thus u satisfies (47) automatically. 
Numerically it is better to work with the change of variables 

w{s) = u(e'*) + 4s, --00 < s < 

which transforms the equation A^u = Ae" into 

iw + 8(A^-2)(iV-4) = Ae™, -00 < s < 
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where 

The boundary conditions u{l) = 0, = then yield 

u;(0)=0, w'{0)^A. 

Regarding the behavior of if as s ^ — oo observe that 

, /8(iV-2)(7V-4)\ ^ 
u(r) — — 41ogr + log I I lor r < tq 



if and only if 



, 8(iV-2)(Ar-4) , „ 

w(sj = log lor all s < logro. 

A 



The steps we perform arc the following: 

1) We fix xo < and using numerical software we follow a branch of solutions to 

' Lw + 8iN -2){N -A) = Xe"^, xq < s < 

w(0) = 0, w'{0)^t 

8{N-2){N-4) d^w <Pw 
w{xo)^\og , -^(xo) = 0, -^(xo) = 

as t increases from to 4. The numerical solution (w, A) we are interested in 
corresponds to the case t — A. The five boundary conditions are due to the fact 
that we are solving a fourth order equation with an unknown parameter A. 

2) Based on w, A we construct a function w which is constant for s < xq 
and piecewise polynomial for ccq < s < 0. More precisely, we first divide the 
interval [xo,0] in smaller intervals of length h. Then we generate a cubic spline 
approximation gfi with fioating point cocfiicients of From g-fi we generate a 
piecewise cubic polynomial gra which uses rational coefficients and wc integrate it 
4 times to obtain w, where the constants of integration are such that ^(xo) = 0, 
1 < j < 3 and w{xo) is a rational approximation of log(8(A^— 2)(A^ — 4)/A). Thus w 
is a piecewise polynomial function that in each interval is of degree 7 with rational 
coefficients, and which is globally C"^. Wc also let A be a rational approximation of 
A. With these choices note that Lw + 8{N — 2)(A^ — 4) — Ae"' is a small constant 
(not necessarily zero) for s < xq. 

3) The condition (44) and (45) we need to check for u are equivalent to the 
following inequalities for w 

(57) Lw + 8(iV-2)(A^-4)-(A + £o)e"' <0, -cxX s < 

(58) iiy + 8(A^-2)(iV-4)-(A-eo)e"' >0, -oo < s < 0. 

Using a program in Maple we verify that w satisfies (57) and (58). This is done 
evaluating a second order Taylor approximation of Lw + 8{N—2){N — 4:) — (A + £)e™ 
at sufficiently close mesh points. All arithmetic computations are done with rational 
numbers, thus obtaining exact results. The exponential function is approximated by 
a Taylor polynomial of degree 14 and the difference with the real value is controlled. 
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More precisely, we write 

f{s) ^Lw + 8{N - 2){N - 4) - (A + eo)e", 

f{s) =Lw + 8{N - 2){N - 4) - (A + eo)T{w), 

where T is the Taylor polynomial of order 14 of the exponential function around 0. 
Applying Taylor's formula to / at yj, we have for s e [yj, yj+h], 

fis)<f{y,) + \.ny,)\h + lMh' 



< hyj) + \~ny,)\h + \Mh' + \f{y,) - ~f{y,)\ + - /'(y,)|/i 

< /(%) + \f'{y3)\h + ^Mh^ + E,+ E^h, 



where 



M is a bound for |/"| in [yj,yj + h] 

El is such that (A + eo)\e"' - T{w)\ < Ei in [yj,yj + h] 

E2 is such that (A + £0)1(6*" - T'{w))w'\ < E2 in [yj,yj + h]. 

So, inequality (57) will be verified on each interval [yj, yj + h] where w is a polyno- 
mial as soon as 

(59) fiy,) + \f'{yj)\h + hlh^ + Ei + < 0. 

When more accuracy is desired, instead of (59) one can verify that 

fix..) + |/'(x.)|- + hii-f + Ei+E-2-<0, 
m 2 m m 

where (xi)i=i...m+i are m + 1 equally spaced points in [yj, yj + h]. 

We obtain exact values for the upper bounds M, Ei, E2 as follows. First note that 

/" = Lw"~{X + eo)e^{{w'f + w"). On we havew{s) = ELo 

and we estimate |if(s)| < X)I=o f'^'' ^ ^ b/jTyj + Similarly, 



(60) 



aw 



< ^ i(i - 1) . . . (i - £ + l)\a,\h^-^ for all s G [yj,yj + h] 



The exponential is estimated by e™ < < 3, since our numerical data satisfies the 
rough bounds —3/2 < w < 1. Using this information and (60) yields a rational 
upper bound M . Ei is estimated using Taylor's formula ; 

(3/2)15 



El = (A + £o)- 



15! 



Similarily, i?2 = (A + £0) ^^(^1 Bi where Bi is the right hand side of (60) when 



14! 

e=i. 



4) We show that the operator — /3e" where u{r) = ?ii(logr) — 41ogr, satisfies 
condition (48) for some P > Po where /3o is given by (49). In dimension iV > 13 
the operator — /?e" has indeed a positive eigenfunction in Hq{B) with finite 
eigenvalue if P is not too large. The reason is that near the origin 
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where c is a number close to 8{N — 2){N — 4)l3/ X. If /3 is not too large compared to A 
then c < N^{N — 4)^/16 and hence, using (6), A^ — /3e" is coercive in HQ{Bra) (this 
holds under even weaker conditions, sec [11]). It follows that there exists a first 
eigenfunction Lpi € Hq{B) for the operator A-^ — /?e" with a finite first eigenvalue 
/^i, that is 

A^cpi — Pe^ifi = Hiifi in B 
ipi> in B 

^1 e Hl{B). 

Moreover /.ti can be characterized as 



Ail = inf ■'^ 



^(A^)2-/3eV 



2 



and is the smallest number for which a positive eigenfunction in H^(Q?) exists. 

Thus to prove that (48) holds it suffices to verify that /xi > and for this it is 
enough to show the existence of a nonnegative tp S Hq{B), f ^ such that 

fAV-/3e"^>0 inB 
ip ^ on dB 



(61) 



^ < on dB. 
on 



Indeed, multiplication of (61) by fi and integration by parts yields 

fii [ ipipi + / ^^^i > 0- 
Jb JdB on 

But Aipi > on dB and thus /ii > 0. To achieve (61) we again change variables 
and define 

0(s) = ip{e'') - oo < s < 0. 
Then we have to find (/> > 0, (/> ^ satisfying 

' - /3e"'0 > in-oo<s<0 

(62) J m = 

^'(0) < 0. 

Regarding the behavior as s — > — oo, we note that w is constant for — oo < s < xq, 
and therefore, if 

L(f> - Pe^'cf) eeO -oo<s<a;o 

then is a linear combination of exponential functions e""" where a must be a 
solution to 

- 2(A^ - 4)a^ + {N^ - ION + 20)a^ + 2{N - 2){N - 4)a = /3e"'^'^°) 
where /3e"'(^«) is close to 8(A^ - 2){N - 4)/3/A. If > 13 the polynomial 

- 2{N ~ 4)a3 + (iV^ _ iQiV + 20)c? + 2(iV - 2)(A^ - 4)a - 8(iV - 2)(N - 4) 

has 4 distinct real roots, while if < 12 there are 2 real roots and 2 complex 
conjugate. If > 13 there is exactly one root in the interval (0, (TV — 4)/2), 2 roots 
greater than (A^ — 4)/2 and one negative. We know that Lp(r) = (pilogr) ^ r~" is 
in H^, which forces a < {N — A)/2. It follows that for s < xq, (p is a. combination of 
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g-aos^ g-Qis ^i^gj-g ao > 0, ai < are the 2 roots smaller than a < {N — i)/2. For 
simplicity, however we will look for cfi such that (/)(s) = Ce^""'' for s < xo where 
C > is a constant. This restriction will mean that we will not be able to impose 
4>'{0) = at the end. This is not a problem because 0'(O) < 0. 

Notice that we only need the inequality in (62) and hence we need to choose 
a e (0, iV - 4/2) such that 

a* ~ 2{N - i)a^ + (iV^ - 107V + 20)a'^ + 2{N - 2){N - 4)a > /?e'"(^"). 

The precise choice we employed in each dimension is in a summary table at the end 
of this section. 

To find a suitable function cj) with the behavior (f)(s) = Ce^" for s < xq we set 
4> = 'ipe~'^^ and solve the following equation 

- I3e^i: = J 

where the operator Tq, is given by 

= ^ + {-4a + 2{N - 4))^ + (6a2 _ 6a(iV - 4) + iV^ - ION + 20)^ 

+ (-4a^ + 6a^(N - 4) - 2a(N^ - ION + 20) - 2(N - 2)(N - 4))^ 

ds 

+ (a^ - 2a^{N - 4) + a^{N^ - ION + 20) + 2a{N - 2){N - 4))?/' 

and / is some smooth function such that / > 0, / ^ 0. Actually we choose P > (5q 
(where /3o is given in (49)) find a satisfying approximately 

~ 2{N - 4)^3 + [N"^ - lOiV + 20)a^ + 2(iV - 2){N - A)a = (ie'"'^''°\ 

We solve numerically 

TaV- - /3e™i/i = / Xo < s < 

^(xo) = 1, i^"{xo) = 0, ij"'{xo) = 

i/)(0) = 

Using the same strategy as in 2) from the numerical approximation of we 
compute a piecewise polynomial -0 of degree 7, which is globally and constant 
for s < Xq- The constant ip{xo) is chosen so that '0(0) = 0- We use then Maple to 
verify the following inequalities 

■0 > xo<s<0 
Tatp - Pe^ip > xo < s < 
^-'(0) < 

where (3q < f3 < j] and < a < {N — 4)/2 are suitably chosen. 

At the URLs: 
http : //www . lamf a . u-picardie . f r/ dupaigne/ 
http : //www. ime .unicamp.br/~msm/ 

we have provided the data of the functions w and defined as piecewise polynomials 
of degree 7 in [xo,0] with rational coefficients for each dimension in 13 < < 31. 
We also give a rational approximation of the constants involved in the corresponding 
problems. 
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Wc use Maple to verify that w and ip (with suitable extensions) are global 
functions and satisfy the corresponding inequalities, using only its capability to 
operate on arbitrary rational numbers. These operations arc exact and are limited 
only by the memory of the computer and clearly slower than floating point opera- 
tions. We chose Maple since it is a widely used software, but the reader can check 
the validity of our results with any other software (see e.g. the open-source solution 
pari/gp). 

The tests were conducted using Maple 9. 
Summary of parameters and results 



N 


A 


eo 




£ 




/3 




a 


13 


2438.6 


1 


5 


10- 


-7 


2550 


2500 


3.9 


14 


2911.2 


1 


3 


10- 


-6 


3100 


3000 


3.4 


15 


3423.8 


1 


3 


10- 


-6 


3600 


3500 


3.1 


16 


3976.4 


1 


1 


10- 


-5 


4100 


4000 


3.0 


17 


4568.8 


1 


2 


10- 


-4 


4800 


4600 


3.0 


18 


5201.1 


2 


2 


IQ- 


-4 


5400 


5300 


2.7 


19 


5873.2 


2 


2 


10- 


-4 


6100 


6000 


2.7 


20 


6585.1 


3 


7 


10- 


-4 


7000 


6800 


2.7 


21 


7336.7 


3 


7 


10 


-4 


7700 


7500 


2.6 


22 


8128.1 


4 


1 


10 


-3 


8600 


8400 


2.6 


23 


8959.1 


4 


1 


10 


-3 


9400 


9200 


2.5 


24 


9829.8 


4 


1 


10 


-3 


10400 


10200 


2.5 


25 


10740.1 


4 


1 


10 


-3 


11400 


11200 


2.5 


26 


11690.1 


6 


2 


10 


-3 


12400 


12200 


2.5 


27 


12679.7 


7 


2 


10 


-3 


13400 


13200 


2.4 


28 


13709.0 


7 


2 


10 


-3 


14500 


14300 


2.4 


29 


14777.8 


7 


2 


10- 


-3 


15400 


15200 


2.4 


30 


15886.2 


8 


2 


10- 


-3 


16600 


16400 


2.4 


31 


17034.3 


10 


2 


10- 


-3 


17600 


17500 


2.3 



Remark 6.4. 1 ) Although we work with A rational, in the table above we prefer to 
display a decimal approximation of X. 

2) In the previous table we selected a "large" value of eq in order to have a fast 
verification with Maple. By requiring more accuracy in the numerical calculations, 
using a smaller value of £o and using more subintervals to verify the inequalities 
in the Maple program it is possible to obtain better estimates of A* . For instance, 
using formulas (50), we obtained 

N ^ £0 £ ^min ^max (3 (3 a 

13 U38.589 0.003 5-10~' 2438.583 2438.595 2550 2510 3.9 

14 2911.194 0.003 5-10-^ 2911.188 2911.200 3100 3000 3.4 

The verification above, however, required to check 1500 subintervals of each of 
the 4500 intervals of length 0.002, which amounts to substantial computer time. 

7. Proof of Proposition 1.5 

Throughout this section, we restrict, as permitted, to the case a = 0. 
a)Let u denote the extremal solution of (1) with homogeneous Dirichlet boundary 
condition a = 6 = 0. We extend u on its maximal interval of existence (0, R). 
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Lemma 7.1. R < oo and u{r) ^ \og{R — r) ^ for r ^ R. 

Proof. The fact that R < oo can be readily deduced from Section 2 of [1]. We 
present an alternative (and more quantitative) argument. We first observe that 



(63) 



i"--u'>0 for all r e [l,i?). 



Integrate indeed (1) over a ball of radius r to conclude that 



(64) 



< A 



9 A 

— Au = ujNr 
or 



N-l I 



N-l , „ 1 

[u u ) 

r r 



If r = 1, since u is nonnegative in (0,1) and u{\) = u'[l) = 0, wc must have 
m"(1) > 0. In fact, u"{l) > 0. Otherwise, we would have u"{l) = and u"'(l) > 
by (64), contradicting m > in (0, 1). So, we may define 

R = sup{r > 1 : u" (t) - ^u'{t) > for aU t e [1, r)} 
and wc just need to prove that R = R. Assume this is not the case, then u"{R) 



{u" - j,u' 



y (R) < 0. This contradicts (64) and we have 



j^u'{R) = and u'"{R) 
just proved (63). In particular, we see that u is convex increasing on (1, i?). 



Since u is radial, (1) reduces to 
.(4) ^ 2(iV-l) . 



(65) 
Multiply by u' 



{N-l)iN-S) „ _ (jV-l)(jV-3) , _ 



r 



5 u u [u ) = A(e j , 



which we rewrite as 

[(u"V)'-u"V'] + 2(7V-l) 



+ {N-l){N-3) 



2r2 



A(e«)' 



By (63), it follows that for r e [1,R), 

'1 



[iu'"u'y - u"'u"] + 2{N - 1) -u"u' + (iV - 1){N - 3) 



2r2 



> A(e")' 



Integrating, we obtain for some constant A 



1 



+ 2(iV- l)-u"u' 
2 r 



{N-l){N-i) [u'f 



> Ae" - A. 



We multiply again by u' 
(66) 



(u" («')')' - (("')')1 - -iu"?u' + 2(N - l)-u"(u')^ 



We deduce from (63) that 



1 



u"{u'f 



1 /I 



2 Vr- 



/\3 



-in') 



iu'n-uy < \ {-{u'f 

r 2 \ r 



and 
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^{ur<kuru"<U-{ur 

j,z J, 2 \r 

Using this information in (66), dropping nonpositivc terms and integrating, we 
obtain for some constant B, 

u"{u')^ + ^^^^—^-(u'f > Xe'^-Au-B 
4 r 

Applying (63) again, it follows that for C = ^4^^ + 1 

Cu"{u'f > Xe^'-Au-B 

which after multiplication by u' and integration provides positive constants c, C 
such that 

(w')* > c(e" - Au^ -Bu- C). 
At this point, we observe that since u is convex and increasing, u converges to +cx) 
as r approaches R. Hence, for r close enough to R and for c > perhaps smaller, 

u' > c e"/\ 

By Gronwall's lemma, R is finite and 

u < —4 log(-R — r) + Cfor rclose to R. 

It remains to prove that u > — 41og(/(! — r) — C. This time, we rewrite (1) as 

We multiply by r^^^(Au)' and obtain : 

i [r2^-2((Au)')']' = Ar2^-2e"(Au)' < Ce"(A«)' < C(e"Au)'. 

Hence, for r close to R and C perhaps larger, 

{{Auyy < Ce"Au 

and so 

VA^(Au)' < Ce^/^Aw < C'e^/^w" < Cie'^/^uJ, 
where we have used (63). Integrate to conclude that 

(Au)3/2 < Ce'^/^u'. 
Solving for Ait and multiplying by {u'Y^^, we obtain in particular that 

Integrating again, it follows that {u')^^'^ < Ce^^^, i.e. 

u' < Ce"/^ 
It then follows easily that (for r close to R) 

u > -41og(i?- r) - C. 

□ 

Proof of Proposition 1.5 a). Given N > 13, let b"^"-^ denote the supremum of 
all parameters b > —4 such that the corresponding extremal solution is singular. 
We first observe that 
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In fact, it follows from Sections 5 and 6 that the extremal solution u associated to 
parameters a = & = is strictly stable : 

67 inf ■'^^ ^' , ^ — — > 0. 

Extend u as before on its maximal interval of existence (0, R). Choosing R £ (1, ^) 
close to 1, we deduce that (67) still holds on the ball Bn. In particular, letting 
v{x) — u(Rx) — u{R) for x Cz B, wc conclude that w is a singular stable solution 
of (1) with a = and b = Ru'{R) > 0. By Proposition 1.7, we conclude that 

f^max ^ -yy^ ^^^^ prOVC that 

5"°^ < oo. 

Assume this is not the case and let u„ denote the (singular) extremal solution 
associated to 6,i, where 5„ y oo. We first observe that there exists p„ G (0, 1) such 
that u'^{pn) = 0. Otherwise, m„ would remain monotone increasing on (0, 1), hence 
bounded above by u„(l) = 0. It would then follow from (1) and elliptic regularity 
that Un is bounded. Let u„(x) = Un(pnx) — u^ipn) for x G B and observe that Vn 
solves (1) with a = b — and some A = A„. Clearly Vn is stable and singular. By 
Proposition 1.7, f„ coincides with u, the extremal solution of (1) with a = b — 0. 
By standard ODE theory, ij„ = u on (0, R). In addition. 



1 

Pn 









(p„) 








Pn 





-OO, 



which can only happen if 1 /p„ R. 

Now, since u„ is stable on _B, m = w„ is stable on Bi/p^. Letting n ^ oo, we 
conclude that u is stable on Bj^. This clearly contradicts Lemma 7.1. 

We have just proved that b™"'^ is finite. It remains to prove that u* is singular 
when —4 < b < 5"*"^. We begin with the case b = b"^""^ . Choose a sequence 
(6„) converging to fe™*^^ and such that the corresponding extremal solution u„ is 
singular. Using the same notation as above, we find a sequence p„ € (0, 1) such 
that 

1 , / 1 
— u 



Pn \ Pn , 

Taking subsequences if necessary and passing to the limit as n — > oo, we obtain for 
some p G (0, 1) 



1 



- =6 



7nax 



— U 

P VP/ 

Furthermore, by construction of u is stable in -Bi/p„ hence in Si/p- This implies 
that V defined for x S i? by v{x) = '"(^) ~ "(p)' ^ stable singular solution of 
(1) with b = Jf^°^^ . By Proposition 1.7, we conclude that the extremal solution is 
singular when b = Jf^"^^ . 

When b = —4 we have already mentioned in the introduction that u* is singular 
for > 13 as a direct consequence of Proposition 1.7 and Rellich's inequality. 

So we are left with the case —4 < 5 < 5™°^. Let denote the extremal 
solution when b = 6™°^, which is singular, and A^ the corresponding parameter. 
For < i? < 1 set 

Ur{x) = U*,{Rx) - <„(i?). 
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Then 

A2wfl = Ai?,e"« where Xr ^ X*QR'^e<^^\ 
and ur = on dB, while 

By (31), note that 

du* du* 

6""^" asi?^l, and R^{R) ^ -A as i? ^ 0. 
dr dr 

Thus, for any —4 < b < 5™°^ we have found a singular stable solution to (1) (with 
a — 0). By Proposition 1.7 the extremal solution to this problem is singular. 

□ 

Proof of Proposition 1.5 b). Let b > —4. Lemma 5.1 applies also for b > —4 
and yields u* < u where u{x) = — 41og|a;|. We now modify slightly the proof of 
Lemma 5.2. Indeed, consider w — {A + 6)(1 — r^)/2 and define u = u — w. Then 

A^u = 8(iV - 2){N - 4)^ = 8(N - 2)(N - 4)e" = 8(iV - 2)(N - 4)e"+™ 

< 8(iV-2)(7V-4)e(4+'')/2e". 

Also u{l) = 0, u'(l) = 5, so u is a subsolution to (1) with parameter Aq = 8(A^ — 
2)(iV-4)e(4+f')/2. 

If is sufficiently large, depending on b, we have Aq < N'^{N — 4)^/16. Then 
by (6) u is a stable subsolution of (1) with A = Aq. As in Lemma 5.2 this implies 
A* < Ao. 

Thus for large enough A^ we have A*e"* < r'^^SiN - 2)(A - 4)e(''+'')/2 < 
r-*N^{N - 4)2/16. This and (6) show that 

which is not possible if u* is bounded. □ 
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